Collisional Granular Flow as a Micropolar Fluid 
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We show that a micropolar fluid model successfuUy describes collisional granular flows on a slope. 
A micropolar fluid is the fluid witlı internal structures in whiclı coupling between the spin of each 
partide and the macroscopic velocity field is taken into account. It is a hydrodynamical framework 
suitable for granular systems wlıiclı consists of particles with macroscopic size. We demonstrate 
that the model equations can quantitatively reproduce the velocity and the angular velocity profiles 
obtained from the numerical simulation of the collisional granular flow on a slope using a simple 
estimate for the parameters in the theory. 
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In spite of the long history of research on the gran- 
ular flows, the theoretical framework for their rheology 
has not yet been established. One factor that makes an- 
alytical treatment difhcult is that there is not a great 
separation of the length scales; the size of each particle is 
often comparable with the scale of the macroscopic col- 
lective motion. Therefore, there are many situations in 
which simple hydrodynamic approaches cannot be used 
to characterize granular flows 0] . Evcn whcn we consider 
the rapid granular flows in which the density is low 
enough that kinetic theory seems to be valid, the coupling 
between the rotation of each particle and macroscopic ve- 
locity field may not be negligible. Thus, the behavior of 
the flow can deviate from that of a Newtonian fluid. 

The micropolar fluid model is a continuum model to 
describe a fluid that consists of particles with spinning 
motion The model equations inciude an asymmetric 
stress tensor and a couple stress tensor. Therefore, the 
model can be a suitable framework to describe granular 
flows. 

Although some research on the application of the mi- 
cropolar fluid model to granular flows has been carried 
out Q, most work has considered dense granular flows, 
in which the constitutivc cquations adopted for the stress 
and the couple stress tensors weıe very complicated. 
Hence it was difficult to interpret the results physically. 

In this paper, we apply a micropolar fluid model to 
a collisional granular flow. We adopt a set of constitu- 
tivc equations that are a simple and natural extension 
of those for a Newtonian fluid. Note that Newtonian 
fluid models of rapid granular flows ||^ have been well 
established. We calculate the velocity and angular ve- 
locity profiles for uniform, steady flow on a slope and 
demonstrate that the micropolar fluid model reproduces 
the results of numerical simulations. 

It is easy to derive the equations for the number density 
n, the velocity «i, and the microrotation field LJi of a 
System that consists of identical particles "vvith mass m 
and moment of inertia /. From the conservation laws of 
mass, momcntum, and angular momentum we obtain 



mnDtVi = mnfi-\-djSj 



InDtOJi 



djCjı 



(2) 
(3) 



The summation convention applies to repeated sub- 
scripts. di represents a partial derivative with respect 
to the ith coordinate, Dt = d/dt + Vk.dk is Lagrange's 
derivative, and fi is the body force per unit mass. Here, 
Sij and Cij are the stress tensor and the couple stress 
tensor that, respectively, represents the j component of 
the surface force and the surface torçue acting on the 
plane perpendicular to the i axis per unit area, and s,-"'' 
is the torque due to the asymmetric part of the stress 
tensor defined as 



(4) 



where e^-fe is the alternating tensor of Levi-Civita. 

For the constitutivc equation of the stress tensor Sij, 
we adopt p|: 



djVi 



- 2eijkUJk] , (5) 



with Sij, Kronccker's delta. The symmetric part of Sij in 
Eq. (|) is the same as the stress tensor of the Newtonian 
fluid with pressurep, shear viscosity fi, and bulk viscosity 
A. The third term on the right hand of Eq. (^) represents 
the asymmetric part of the stress tensor due to the differ- 
ence between the rotation of the mean velocity field and 
particles' own spin. This gives s'-"-' = 2/ir[V x t; — 2lü], 
using Eq. . The microrotation viscosity represents 
the coupling between the velocity and the microrotation 
field. For the couple stress tensor dj, we use the theorem 
that isotropic second order tensors can be expanded as 
the trace, the symmetric part, and the asymmetric part 
of rate of strain tensor for microrotation g : 



Cij = codkUJkStj + 



{diLüj + djUJi) 



(6) 



Dtn = -ndkVk, 



(1) 



where the coefîicients of angular viscosity are cg, fJ.A, and 
fiB- It should be noted that the dimensions of these 
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coefHcients are different from the viscosities in the stress 
tensor by length-squared because of the difFerence of the 
dimensions between Sij and Cy and between Vi and uJi. 

It is debatable whether such a straightforward exten- 
sion of Newtonian constitution relations can be apphed to 
granular flow, because the hydrostatic term in a granular 
material should have a different form. For a coUisional 
flow, however, we expect that such an effect is not impor- 
tant. Therefore, we concentrate on a colhsional granular 
flow. 

The coefîicients of viscosity that appear in Eqs. (||) 
and (|^) have been derived based on the kinetic theory 
of three-dimensional spheres with rough surfaces jşl I). 
Here, for later convenience, and also to make the phys- 
ical meaning of the model clear, we briefly summarize 
the rough estimate of the coefficicnts of viscosity to the 
lowest order for two-dimensional disks using elementary 
kinetic theory. Let us consider a two-dimensional fluid 
that consists of idcntical disks with diameter a and that 
is flowing uniformly in the x direction, namely n — n(y), 
V = (u(y),0, 0), and lj = {0,0,uj{y)). Thcn we have 



Syx = Mw'(y) + fJ-r[u'{y) + 2[j(2/)], 



and 



fJ'Bi^'iy), 



(7) 



(8) 



where the prime denotes a derivative with respect to its 
argument. Here, Syx (Cyz) is the x (z) component of the 
force (torque) per unit area acting on the plane perpen- 
dicular to the y axis. 

The coefficient /i in Eq. (|^) corresponds to the kinetic 
viscosity in dilute gas, which we can find an estimate 
of kinetic theory in textbooks on statistical physics, e.g. 
Ref. M. It is given by 



(9) 



where l is the mean free path and v is the mean square 
displacement of the velocity in the y direction. Here, T 
is the granular tcmperature and the relation l ~ l/{na) 
in two-dimensional system is used. 

The coefficient /is is estimated by a similar argument 
to that of /i. Because Cyz represents the net angular 
momentum transfer per unit time per unit length due 
to the microrotation by particles crossing the plane y = 
const, we can use the argument for ^ by replacing u{y) 
and m by and /, respectively. Then /zb is estimated 
as 



MS 



nvll 



(10) 



with / = TO(T^/8 for a two-dimensional disk. As we have 
mentioned, the dimensions of /z and fiB are different. 

For an estimation of ^r, which gives coupling between 
particles' own spin and velocity field, we consider the 



collision of two disks near the plane y = const . If the 
surface of the disk has some roughness, the momentum 
tangent to the relative position of the coUiding particles 
at the time of contact is transferred from one partide 
to another. It is plausible to assume that the tangential 
momentum transferred in one collision is proportional 
to m times Au, the relative tangential velocity of each 
partide at the contact point. In order to simplify the 
estimation, let us consider the situation with a uniform 
velocity field, namely u'{y) = O Then Au is given by 
Au ^ 2{a /2)uj{y). Because the frequency of collision per 
unit time per unit length near the plane y = const with 
the width a is proportional to n^ücr^ in two-dimensions, 
the momentum transfer across the plane by coUisions is 
estimated as 

AM r~.n^va^m[2uj{y)]. (11) 

Comparing Eq. ( |IÎ]) and the second term of Eq. (0) with 
iJ-'iy) = Ol obtain 

/ir ~ n^va^m ~ n^a'^VTm. (12) 

Summarizing the results above, which are consistent 
with Ref. wc gct following expressions for the cocffi- 
cients of viscosity; 



/ir 



(ncr2)2 



(13) 



It should be noted that, because the dimcnsion of /is is 
different from that of /z and /i^, we need to introduce 
another length scale to characterize the macroscopic flow 
behavior in order to compare them. On the other hand, 
we can see that, when the number density is high enough, 
/il- becomes comparable to /i, then the coupling between 
the angular momentum and the linear momentum should 
play an important role. 

Now we present the uniform, steady solution of the mi- 
cropolar fluid equations (|î|), (||), and (^ on a slope, and 
compare the obtained profiles with the result of numcri- 
cal simulation. 

Let us consider the two-dimensional steady flow on 
a slope under the gravity. We take the x axis in the 
direction down the slope, and the y axis in the direc- 
tion perpendicular to the slope. The inclination an- 
gle of the slope is 6, and the acceleration of gravity 
g — {g sm9, —g cos9,0). In the uniform, steady flow, 
mn ^ p ^ p{y), v = (u(y),0,0), and a; = (0,0,w(y)), 
therefore the equation of continuity (|î|) is automatically 
satisfied. Eqs. (||) and (^) reduce to following diffcrcntial 
equations: 



pgsmO+ — 
dy 



du 
^ dt/ 



dit 

-- + 2üj 
dy 



-pgcosö 



dp 
dy 



= O, (14) 
= O, (15) 
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FIG. 1: The properties of the uniform flow. FiUed circles 
show the data from the simulation and the sohd hne shows the 
uniform, steady solution of the micropolar fluid equations. (a) 
The angular velocity profile Lü{y). The dashed hne is the finite 
series approximation of the solution Eq. (jİ^) with fc = 1 to 
5; the difference from the solid line is hardly distinguished by 
the eye. (b) The profile of the deviation of the microrotation 
from the rotation of velocity field, x v — uj. 



Its general solution is given as the sum of a particular 
solution üjp and two homogeneous solutions lui and â'2 
by 



Û{Y) = ACj^İY) + BCJ2İY) + Cjp{Y), 



(22) 



with integration constants A and B. Changing the vari- 
able from y to 77 = exp(— y), we can obtain the expres- 
sions 



Lüp 
a.2k+l 

b2k 



k=l 

f2k+ia.2k-l, 



2k+l 



with 03 = 



2a 



Y^b2kV^^, with 60 = 1, 

k=0 



(23) 



(24) 



and 



and 



2/ir 





d 


da> 




dy 


dy 



= 0. 



(16) 



The equation of state for a two-dimensional granular gas 
has been derived in Ref. jî^, but here we adopt the 
lowest order estimate, namely p = pT/m. With the aid 
of the assumption of constant temperature T = T, we 
obtain the density profile from Eq. (^5|); 



p = po exp 



(-!) 



T 



mg cos 9 



(17) 



Then, from the estimate of the coefficients Eq. (13), we 



introduce the non-dimensional constants a and f3 as 

aexp(-2î;//ı), ^ = /îcr^. (18) 



After integrating Eqs. ( [î^ ) and (^ with Eqs. ( p7| ) and 
( p^ ) , we obtain the following relation between the veloc- 
ity field u{y) and the microrotation field Lü{y)] 



Ü{Y) = - exp(-y) 



dCj{Y) 

~2 dr 



where 



Pogh^ sinö//x' 



Pogh sind / 



+ Ao, (19) 



(20) 



with Y = y/h and e — a/h. Here, ^0 is an integration 
constant. In the integration of Eq. (|ÎJ), we imposed the 
boundary condition for stress tensor at the free surface, 
limj^^oo Syx 0. From Eqs. (16) and (p^, we have 



d^Cj{Y) ^aexp{^2Y) 



dF2 



exp(-r) -f 2Ü{Y) 
1 + Q;exp(-2r) 



= 0. (21) 



= â)ı logT; + ^ C2fc77^'', with co = O 

fc=l 

a r 1 

fc3 



UI2 



C2k = /2feC2fe-2 



k-1 



where 



fn 
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2 - (m- 2) 



b2k-2, (25) 



(26) 



Note that ûp and (Di are regular but 0)2 diverges at 77 = O 
{Y 00). 

Now we show that the solution obtained from the mi- 
cropolar fluid model can quantitatively reproduce the re- 
sults of numerical simulations of coUisional granular fLow 
on a bumpy slope ||îî|. In the simulation, the discrete 
element method is employed with normal and tangential 
elastic force and dissipation. The particles are modeled 
by disks with diameter ct, and the parameters are cho- 
sen so that the normal restitution e„ becomes 0.7. The 
surface of the slope is made rough by attaching identical 
particles to it. In the following, ali quantities are given 
in dimensionless form with mass unit m, length unit a, 
and time unit r = \föjg. 

We used the data of the simulation with inclination 
sinö = 0.45, system size L = 1002 and number of fLow- 
ing partide N = 1000. In the simulations, uniform flow 
is realized during 500 ^ t ^ 2000 where t is the time . 
In order to determine the profiles of the mean quantities 
describing the flow, we divide the space into layers which 
are one partide diameter wide parallel to the slope, calcu- 
late the averages inside layers, and then average the data 
över the time within the uniform flow, 1000 < t < 1500. 
The origin y = O is taken one diameter above from the 
top of the disks attached to the slope. 

Now let us compare our analytical result from the mi- 
cropolar fluid model with the simulation data. First, we 
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confirmed that the number density profile can be well- 
fitted by Eq. ^) with h = 2.24, namely f - 2.0, which 
has been checked to be close to the averaged value of T. 
In order to fit the solutions and (^2|) to the simu- 
lation data, we determine the boundary values u(0) and 
uı{0) from the data and treat a, (3, fi/po and uj'{0) as 
fitting parameters. From Fig. ^ (a), we can see that 
the micropolar fluid equations reproduce the angular ve- 
locity quantitatively. Here, the value of parameters are 
a ~ 0.10, /3 ~ 0.12, /x/po ~ 0.95, and w'(0) ~ 2.9. The 
velocity profile u{y) can also be well-reproduced by the 
solution (^; however, because the density is low, the 
deviation of the velocity profile from the solution of the 
Navier-Stokes equation, i. e. ü{Y) = — exp(— F) + ^o, is 
small. Actually, we have checked that the data can be 
reasonably fit with any small value of a, as long as a/P is 
chosen appropriately. However, it is significant that the 
solution (|2^ ) can reproduce the sharp variation of üj{y) 
över only a few diameters near the base. 

The more important thing is that the mean spin u: 
deviates systematically from the rate of bulk rotation 
İV X ı>, and that the micropolar fiuid model can re- 
produce this deviation. In Fig. (b), we see that the 
deviation is large near the base, because each partide 
is forced to rotate by the coUision with the slope. This 
result indicates the importance of the couple stress near 
the boundary, as pointed out in Ref . |jî^ . This deviation 
may produce the velocity profile different from the New- 
tonian fluids described by the Navier-Stokes equation due 
to the coupling between the spin of each partide and the 
linear velocity field On the other hand, it seems the 
deviation also becomes large in the region far from the 
boundary. The reason is that the microrotation field in 
this region is dominated by a small number of particles 
spins that are generated by coUisions with the boundary. 

In summary, the micropolar fluid model has been ap- 
plied to a coUisional granular flow with relatively low 
density. It has been demonstrated that the solution for 
uniform, steady flow on an inclincd surface reproduces 
the results of numerical simulation. 

Because the density is low, the ideal gas assumption 
for the equation of state and the estimation of the vis- 
cosity from the elementary kinetic theory work well. In 
order to apply the model to denser coUisional fiows, we 
need a systematic extension of the theory using, for ex- 
ample, Enskog theory. This has already been done in the 
context of polyatomic fluids for completely rough spheres 
without energy dissipation |^ and it should be possible 
to extend such results to the dissipative case The 
equation of state for dense granular gas has also been 
discussed recently [0. On the other hand, it is known 
that the effect of the partide spin becomes more impor- 
tant in dense frictional flows ||î^. Much research has 
been done to construct the mechanics of granular niedia 
which is valid not only for coUisional fiow but also for 
denser situations based on the micropolar or Cosserat 



theory jl^]. The concepts of micropolar mechanics may 
help to link the understandings of the coUisional and the 
frictional flows of granular materials. 
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